General properties of the pion production reaction in nuclear matter 
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dependence of the reaction, which was measured using 2 H, 4 He, 12 C, ie O, 40 Ca, and 20S Pb targets 
at ~280 MeV. Some general properties of the pion-induced pion production reaction in nuclear 
matter will be presented, based on the combined results of the two studies. 
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1. Introduction 

The 7r + — > ty+'k^ reactions {n2ir) on 45 Sc were studied at intermediate energies, T n + 
= 240, 260, 280, 300, and 320 MeV. The A-dependence of the ix2ix process was studied 
previously on 2 H, 4 He, 12 C , 16 0, 40 Ca and 208 Pb nuclei at T w+ -280 MeV [H El El El U\ ■ 
The results of those investigations, combined with the present investigation of the tt2tt T- 
dependence, allow for the determination of some general properties of the tt2tt reaction in 
nuclei, as well as those of the titt system. The 2 if studies were performed both to understand 
the 7r27r behaviour on a neutron and a proton by means of the ttN — > ttttN quasifree reaction, 
and to observe medium modifications of the mr interaction through the direct comparison to 
the 7r27r data in nuclei. To ensure a reliable comparison, all the 7r27r data were taken under 
the same kinematical conditions. 

The detection of titt pairs in the 1=0 J=0 channel is a way of examining the existence of 
the light scalar-isoscalar meson, the a meson. The direct observation of cr's in vacuum has 
always been frustrated by the elusiveness of the meson, perhaps due to the large width of its 
spectral function jHj. The threshold behaviour of the ( y 7m)i = j = o system has been the subject 
of experimental investigations: the invariant mass intensities display a remarkable strength 
in nuclei, while they appear depleted in the nucleon (i.e., in vacuum) jjj. An interpretation 
of this observation relies on the partial restoration of chiral symmetry in nuclear matter. In 
a chiral symmetric space, chiral partners would have same spin but opposite parity, and be 
degenerate in mass. In this framework, the chiral partner of the a (J p = + ) is the pion 
(J p = _ ). These mesons, however, have the masses which are a few hundred MeV apart. 
The shift of the tttt invariant mass toward the lower 2m n in nuclei is taken as a signature of 
partial restoration of the broken symmetry [HIGH]. ^ ^ s worthwhile noting that this process 
is observed in nuclei, i.e. at nuclear densities p < p n the saturation density. In the same 
environment, another process contributes to the reshaping of the (7r7r)/ = j =0 spectral function 
at around the 2/71^ threshold (and below it), which is due to the P-wave coupling of pions 
to particle-hole (p-h) and delta-hole (A-/i) configurations. The net result is similar, although 
less pronounced, to that of partial restoration of chiral symmetry [TTJ H21 HH1 UH UH] • 

Correlated pion pairs have also been studied in a unitary chiral nonperturbative model 
[16 . In this approach, the a resonance in vacuum is dynamically generated by the strong tttt 
interaction. In nuclear matter, the a properties are only weakly modified ^7j, whereas the 
(tt7t)j = j = o interaction is altered by the P-wave coupling of pions to p-h and A-h states. In 
the 1=0 J=0 channel, ImT n7T enhances its strength at 2-171^ as pion pairs probe increasing p's. 
The increase is detectable at p < p n . The pion production reaction in nuclei proceeds via 
the quasifree nN — > itttN elementary process JHj, and the nuclear influence on the process 
is accounted for by a complete model |15j . which is considered the state-of-the-art in this 
field. 
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Table 1: Characteristics of the particle beam: Central momentum (p, [MeV/c]) and energy 
(T, [MeV]), momentum spread (5p/p at FWHM, [%]), intensity (N, [10 6 particles/s]), and 
composition (vr[%],p[%]). 



P, T 


5p/p 


N 


7T, p 


353, 240 


2.0 


1.6 


95.5, 4.5 


374, 260 


1.1 


3.4 


94.3, 5.7 


396, 280 


1.1 


2.6 


94.9, 5.1 


417, 300 


1.0 


2.4 


95.7, 4.3 


438, 320 


1.6 


1.5 


96.4, 3.6 



The purpose of the present article is to give an account of new results on the T-dependence 
of the 7i2n reaction, and to examine the data in relation to the above mentioned theoretical 
results. In addition, complemented by the earlier studies of the A-dependence, general 
properties of the n2n process in nuclear matter are sketched out. The n2n data were collected 
under the same kinematical conditions to allow for a straightforward comparison of the 
results. Pion pairs were analysed down to opening angles ~ 0° in order to determine the nn 
invariant mass at the 2171^ threshold, and the energy of incident pions was limited to 320 
MeV to reduce the final state interactions of pions with the residual nucleus. The article 
is organised as follows: some features of the experiment are reported in Sec. 2. Sec. 3 
deals with the data analysis. Exclusive tt2tt data at 300 MeV are discussed in Sec. 4. The 
7r27r model used for comparison is presented in Sec. 5. The T-dependence of the reaction is 
reported in Sec. 6. Other existing results are mentioned in Sec. 7. Finally, conclusions are 
summarised in Sec. 8. 

2. The experiment at TRIUMF 

The aim of the experiment was the measurement of differential cross sections for the 
7r + — > 7i + 7i ± reactions at several incident pion kinetic energies. This required single particle 
mass-identification and vector momentum determination for each event. 

2.1 The pion beam and the beam counting 

The experiment was performed at the TRIUMF Meson Facility, using positive pions from 
the Mil channel. The final focus was located at the centre of the detector. Pion energies 
were 240, 260, 280, 300, and 320 MeV. Positive pions were accompanied by the 7r + 's decay 
products, i.e. /x + 's and e + 'e, and p's. However, the contamination of positrons and positive 
muons was found to be less than 1% at all energies. Protons were initially separated from 
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pions by means of an absorber placed in the mid-plane of Mil, and finally intercepted by 
the slits of the channel. Protons escaping the slits were constantly monitored through their 
time-of-flight along the beam line (~13.5 m long). The beam characteristics are summarised 
in Table 1. The target used was a solid slab of 45 Sc 0.558 g/cm 2 thick. The target area 
exposed to the pion beam was 4.25x3.03 cm 2 , which could contain a large fraction of the 
beam envelope. 

The number of pions impinging upon the target was constantly monitored by an in-beam 
scintillator counter, which was located at about 0.8 m from the target. Therefore, the real 
number of pions arriving at the target was corrected for the rate of pion decay, which was 
~3.5% at ^+=280 MeV. In addition, a rate correction due to multiple pions per beam 
burst was applied; this was, for example, ~4.9% at ^+=280 MeV. Detailed calculations are 
reported in Ref. [7j. 

2.2 The CHAOS spectrometer 

Multiparticle charged events were detected by the CHAOS spectrometer [TT?]. CHAOS 
consists of a dipole magnet, four cylindrical wire chambers (WC) and a cylindrical 



Figure 1: Reconstructed particle trajectories in CHAOS for irf — > tt+tt p at T^+^300 MeV. The experi- 
mental layout shows the geometrical disposition of the wire chambers (WC), the first level trigger hardware 
(CFT) and the magnet return yokes in the corners. The CFT segments which are hit by particles are marked 
with crosses, and the energy deposited in the first two layers (AE1 and AE2) is indicated by boxes. AE1 
is hit by two pions, however the thinner segmentation of AE2 is sufficient to distinguish them. 
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telescope of fast counters (CFT). Fig. 1 shows the geometrical disposition of the magnet 
return yokes, wire chambers and telescope, and the trajectories of a reconstructed event, 
71"/" — > 7T + 7T~p on 45 Sc at T 7r +=300 MeV. The magnetic field was varied with T n + to maintain 
constant the trajectory of the incident beam. The field was 0.55 T at T 7r +=240 MeV and 
0.69 at T 7r +=320 MeV. The inner and the outer WC's have a diameter of ~23 cm and ~123 
cm, respectively. The latter is a vector WC, which operates in the fringe field of the magnet. 
The overall WC efficiency to fully reconstructed 7r27r events about 50±5%. With this setup, 
the vector momentum p* n (p, 0) of a 130 MeV/c pion was analysed with 5p~ 5 MeV/c and 
50 < 2°. Both uncertainties are due mainly to pion multiple scattering. The CFT telescope 
encircles WC4. The CFT hardware is composed of two layers of NE110 plastic scintillator, 
0.3 cm and 1.2 cm thick respectively, and a layer of SF5 lead-glass of ~5 radiation lengths. 
The telescope is segmented. Each segment covers an azimuthal angle A© = 18°. Segments 
can be removed to allow the pion beam to enter and exit CHAOS. The first CFT layer 
subtends the smallest zenith angle A$=±7°; therefore, it defines the geometrical solid angle 
of CHAOS Q = 1.5 sr. The necessity of removing 2 (out of 20) CFT segments during 
a measurement leads to a nonuniform acceptance of CHAOS. In addition, pions decaying 
inside CHAOS further increase the irregularity of the acceptance. Thus, GEANT Monte 
Carlo simulations were required, which resulted in assigning a weight to each tc2tt event 
before it was binned. The weight distributions resemble those reported in Ref. [IJ Fig. 4: 
they peak at weight~1.5 and monotonically decrease to zero at weights . In order to avoid 
large corrections in the distributions, the soft cut < weight < a + 2/i on reconstructed 
events is applied, where a and [i are the mean value and the standard deviation of the weight 
distributions, respectively. In the case of T„-=300 MeV, cr + 2/x ~ 5.5. 

2.3 Particle mass-identification 

Incident pions produce background reactions whose rate overwhelms the tt2it rate. There- 
fore, to study the tt2ti process in nuclei, a thorough identification of particle masses is re- 
quired. In the case of CHAOS, this is accomplished by binding the momentum (polarity) of 
a particle to the pulse- height response of the telescope layers. In the present measurement, 
pion momenta do not exceed 230 MeV/c and the particles involved are 7r's, e's, p's and ofs. 

The pulse height response (PH ) of pions and protons in plastic scintillators is PH P / PH V > 
5 for p,r < 230 MeV/c see Fig. 8 in Ref. [2U] and Fig. 3 in Ref. [3]. By using only the first two 
CFT layers, pions were separated from protons (deuterons) with high (~100%) selectivity. 
The same technique cannot favourably be applied to 7r's and e's for momenta exceeding 110 
MeV/c. Thus, SF5 lead-glass Cherenkov counters were used. 

Electrons (positrons) are byproducts of single charge exchange processes in the 45 Sc 
target: 7T + — > 7r° followed by 7r° decay 7r° — ► 77 and 7 conversion 7 — > e + e~. Thus, e + e~ 
pairs may look like 7r + 7r~ pairs. Fig. 2 shows the opening angle distributions of e + e~ (O e + e -, 
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full squares) and 7r + 7r~ (Q w + n - , full diamonds) pairs for the 300 MeV run, where G is defined 
in the lab frame as the difference of the azimuthal angles of reconstructed pairs. Most of 
the e + e - strength occurs near 0°. An e + e~ to 7r + 7r~ misidentification would result in a 
contamination of Q n + n - at around 0°, and thus in a contamination of the tt + 7t~ invariant 
mass at around the 2m n threshold (Ref. [Zj Fig. 13). In order to safely separate electrons 
from pions, the following strategy was followed. Soft kinematical cuts were assigned to reject 
those 7r27r spurious events whose momenta and total energy sum exceed the values allowed 
by the A5 Sc(tt + , tt + tt^p) 4:A Sc phase space. The same soft cuts were also applied to e + e~ 
pairs, which brings about the distributions in Fig. 2. The Cherenkov counters have a 7r to e 
discrimination efficiency above 95% for momenta below 230 MeV/c (Ref. [20] Fig. 11), which 
brings the e + e~ to 7r + 7r~ misidentification rate to about 0.25%. This is further reduced to 
about 0.1% by rejecting events with opening angles less than 3°. For the distributions shown 
in Fig. 2, there are about 3000 e + e - events between 0° and 18°, which is reduced to about 
3 events by the CFT filtering combined with the 3° soft cut. These events would mainly 
affect the 5° 9^+^- datum which consists of about 900 events. Thus the surviving e + e~ 
contaminants have a negligible impact on the opening angle distributions, and in general 
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Figure 2: Opening angle distributions of e + e~ (full squares) and tt + it~ (full diamonds) pairs at T ff +=300 
MeV. The inset diagram shows the tt + tt~ opening angle distribution (open diamonds) for the 7r + — ► tt + tt~p 
reaction; in this case, tt + tt^ pairs are detected in coincidence with protons. 



on the 7r27r spectra. This feature is common to all the examined energies. In the upper 
panel, of Fig. 2, the 0^+^- distribution (open diamonds) is shown for pion pairs detected 
in coincidence with protons, namely, for pairs from the tt + — > tt + tt~p reaction. In this case, 
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kinematical cuts on the kinetic energy sum of particles detected in the reaction exit channel 
(Tsum) avoid contributions from e + e~ pairs in coincidence with a proton. Nevertheless, the 
two ©tt+tt- distributions show a similar behaviour. This second self-reliant result confirms 
that e + e~ events are efficiently rejected. 



3. Analysis 



The 7r27r data are fully corrected within the geometrical acceptance of CHAOS. No at- 
tempt to correct for data outside the CHAOS acceptance has been made, since this would 
involve relying either on phase space simulations or on model calculations. Such a model- 
dependent approach would alter the shape of distributions; in fact, the data must be extrapo- 
lated over a large $— interval without the assistance of experimental data. This is illustrated 
in Fig. 3 for the ittt invariant mass distributions at T 7r =280 MeV: the full-line diagrams are 
the results of phase space simulations for the A5 Sc(n + , 7r + 7r ± A^) 44 X reactions for an ideal 
An detector, while the shaded diagrams depict the effects on distributions of the CHAOS 
geometry. In the present work, the shape of distributions will reflect the coplanar geometry 
of the spectrometer. 



7T -> 7T 7T 




270 310 350 390 430 270 310 350 390 430 

M (MeV) 



Figure 3: Results of phase space simulations for the 45 Sc(n + , 7r+7r ± N) 44 X reactions at 280 MeV. Full- line 
diagrams: tttt invariant mass distributions for an ideal 4tt detector; shaded diagrams: ttt: invariant mass 
distributions when including the CHAOS geometry. The diagrams are normalised at low and high invariant 
masses where distributions appear to have the same behaviour. 



Once the particle mass is identified, the particle vector momenta are checked against 
simulations of the 45 5c(7r + , 7r + 7r ± A^) 44 X phase space, which include the CHAOS acceptance 
and resolution. A reconstructed event is rejected when it exceeds the allowed phase space 
volume. This procedure is applied to all the five energies. As an example, for the 300 MeV 
run the following soft cuts are applied: p n + <215 MeV/c, p n - <215 MeV/c and Tsum <156 
MeV. As earlier mentioned, each ir2ir event is weighted before being binned. In order to 
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avoid large corrections, the soft cut < weight < 5.5 on reconstructed events is applied. 
In the case of the 7r + — > 7r + 7r~ channel, the number of reconstructed events which passed 
the above tests are several thousands for each energy. The number drops by an order of 
magnitude for the 7r + — > 7r + 7r + channel. 




30 60 90 120 150 180 30 60 90 120 150 180 
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Figure 4: Out-of-plane dependence of da/d^^-ir as a function of $„■„■ = $7r lW2 = — for 

reaction channels at T 7r +=280 MeV. Full diamonds: experimental points; full lines: 



the 



r+,r± 



curves linearly joining the experimental data; full histograms: results of phase-space simulations for the 
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N) 4A X reactions normalised to the data at 0° ± 7° and 180° ± 7°. 



The data reduction is based on n — > tx\tx2 events which passed the following tests: have 
the origin in the target region, be mass- and charge-identified, survive the tests on kinematics 
and weights. Many-fold differential cross sections d a / (dTdQ) 7ri (dTdQ) 7T2 are then formed, 
where T is the kinetic energy of a final pion and Q is the solid angle into which a pion 
is scattered. Q n is related to the zenithal angle Q n and the azimuthal angle via the 
equation d£l n = dcos(Q n )d& n . The two angles reflect the coplanar geometry of CHAOS, 
which restricts the out-of-plane acceptance to vary within ±7°, while particles are accepted 
from 0° to 360° in the plane of the reaction except for two angular segments each 18° 
wide. The latter account for the removal of two CFT's which allow the pion beam to enter 
and exit CHAOS (see Fig. 1). The data are represented as triple differential cross sections 
d^a/dO^dQ^dQ^, where O n7T represents (T or Q)^ or a combination of them, and 7r denotes 
a charged pion. d Q^ dn is determined by the quantities f e N ^ , where f e is a parameter which 
is determined by the experimental conditions, N(0) is the number of weighted events in a 
given bin and AO is the bin width for the O observable. The systematic uncertainty in 
assessing f e is 13.4% (er) at the central energy T 7r +=280 meV. A more detailed analysis for 
T„.+ =280 MeV is reported in Ref. [7j. The error bars reported on the ir2iT spectra account 
only for the statistical uncertainty. The overall uncertainty is obtained by summing the 
systematic and the statistical uncertainties in quadrature. In order to obtain the total cross 
section ot, the differential cross section d z a / dO^dVL^dVL^ is summed over and fl^s. 
Fig. 4 illustrates the out-of-plane dependence of the cross section dcr/d& n7T at 280 MeV, 
where <& wn = $^ 17ra = 1$^ — $^1- The integration over the unmeasured portion of & n7r 
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is accomplished either by using a linear function joining & W7T = 0° to ^ n7r = 180°, or by 
interpolating the two measured points with the reaction phase space. Such an approach is 
applied to all the incident pion energies and the results are similar. The values of <jt so 
obtained are listed in table 2 (Sec. 6.4). 



4. The 7T + — > 7r + 7r p exclusive measurement at T 7r +=260 and 300 MeV. 



At 7^+ =300 MeV kinetic energy, the txItx data were collected with sufficient statistics 
(over 10 4 reconstructed events) to enable an analysis of the n + {n) — > it + it~p exclusive reac- 
tion channel. Some results are discussed below. 



Fig. 5 shows the missing energy (E M ) distribution of the ir + 45 'Sc — > it + it~p 44 Sc reaction 
at T 7r += 260 and 300 MeV. The missing energy is defined by the equation Em = [M^sc — 
Muse — m n — mp] — [Tl+ + + T£ — T*+], where the symbol M{m) indicates a nucleus 



7T 



7T 7T p 



T -300 MeV 




T =260 MeV 




-100 -50 50 100 

Missing Energy (MeV) 



7T+7T p 44 Sc reaction at incident kinetic energies 



Figure 5: Missing energy distributions of the 7r + 45 Sc 
of 260 and 300 MeV. The definition of missing energy is given in the text; in this case, Tr + ir~ pairs are 
detected in coincidence with protons. The curves are a Gaussian fit to the data. 



(particle) mass, while % and / denote the initial and final particles, respectively. In this 
framework, Em represents the excitation energy of the 44 Sc residual system. The Em mean 
values, obtained with a Gaussian fit to the data, are ~7 MeV at 7^+ =300 MeV, and ~3 
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MeV at T 7r +=260 MeV. These values are consistent with Em ~0 MeV, when considering 
that the residual nucleus may be left in an excited state. Therefore, the 7r — > -kit reaction at 
intermediate energies is a quasifree process, which mainly involves a single nucleon ttN — > 
7T7tN. A similar result was found in Ref . [3 j , where the n2n reaction was studied as a function 
of the nuclear mass number. 

Fig. 6 shows the invariant mass distributions (shaded diagrams) of tt + tt ± pairs that are 
detected without the requirement of a proton in coincidence. The M niT distributions are also 
plotted as a function of p wn = \p n + +p n ±\, in order to illustrate the p^— dependence of 
for the two reaction channels. The two channels have a similar behaviour: the low M W7T bins 



Jp nj > L75 MeV/c 
£ P™ £ W3 MeV/e 









■ 


" 




;;;;!. 
■ * 



250Q ■ 

j £000 ■ 

1500 ■ 
1000 



£ p m £ L7B MeV/c 



M {MeV) 



Figure 6: Invariant mass distributions of tt + tt ± pairs at 300 MeV for 45 Sc as a function of p^ir, the 
magnitude of the pion pair momentum. Shaded diagrams < p ff7r < 350 MeV/c (the maximum momentum 
reached by a pion pair), full squares < p V7r < 175 MeV/c, and full diamonds 175 < < 350 MeV/c. 

are mainly populated by pairs (full diamonds) with the highest momenta p nn >175 MeV/c. 
These pairs, however, show a different behaviour when their opening angle distributions are 
compared. Fig. 7 depicts the Q nn distributions for p nn >175 MeV/c and 2m n < M nn < 310 
MeV. In this invariant mass range, 7r + 7r + (open diamonds) and 7r + 7r _ (full diamonds) pairs 
are known to carry 1=2 J=0 and 1=0 J=0 quantum numbers, respectively The 
9 vr + 7r + distribution is well described by the n + Ab Sc — ► -n + -n + n 4A Ca phase space (full line). On 
the other hand, the 0^+^- distribution displays a sharp peak at around 0°, before continuing 
as 7r + 45 Sc — > 7T + 7i~p 4A Sc phase space (double line). For both channels, simulations include 
the CHAOS acceptance as well as the p nn and M n7T cuts, and are normalised to the data from 
©tttt > 25°. In the case examined, the only difference between 7r + 7r + and 7r + 7r~ pair is their 
quantum numbers, demonstrating the strong influence of nuclear matter on the (7r7r)/ = j =0 
interaction. Low M n + n - are populated by (7r7r)/ = j =0 pairs which decay preferentially with 
small 0,^ and high p wn . 
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Figure 7: Opening angle distributions of tt + tt + and 7r + 7r _ pairs at 300 MeV, for 175 < p^ < 350 MeV/c 
and 2m„ < < 310. The full and double lines are the results of phase space simulations for the reactions 
7r + 45 Sc — ► Tr + n + n iA Ca and 7r + 45 Sc —> ir + ir~p 44 Sc, respectively. Simulations include the CHAOS 
acceptance as well as the p n7r and M n7r cuts, and are normalised to the data from <d n7r > 25°. 

5. Model of the n2ix reaction in nuclear matter 

The experimental results will be compared with the model predictions of Ref. ^3]. There 
are several good reasons for doing so. The calculation of the n2ix cross sections is based on 
a Monte Carlo method, which permits the CHAOS acceptance to be easily implemented in 
the code. The code requires the energy of the incident pion as an input parameter, thus the 
experimental distributions can be compared with the cross sections precisely at the same 
energies. For each accepted event, the code returns the vector momentum of each pion. 
A similar approach is followed in the data analysis, facilitating the comparison between 
calculated cross sections and distributions of experimental observables. 

The model describes the elementary pion production reaction by involving only one nu- 
cleon of the nucleus, via the quasifree process tcN — > 7r7riV jT3 CHI- The construction of 
the ttN — > ttttN amplitude relies on one-point diagrams, which include both the pion-pole 
and contact term, and two- and three-point diagrams. Nucleons as well as the A(1232) and 
iV*(1440) can be excited as intermediate states. Standard nuclear effects, which may modify 
the reaction cross section, are accounted for. These include Fermi motion, Pauli blocking, 
pion absorption and quasielastic scattering. The role played by the nuclear medium on 
(vT7r)/ = j =0 interacting pairs is carefully examined. A tangible modification of the 7T7T inter- 
action is determined by the coupling of pions to p-h and A-h states, which enhances the 
intensity of the tttt spectral function at around the 2m n threshold. The effects of the nuclear 
medium on the a meson properties are analysed in separate articles fOJEI- The a mesons 
are dynamically generated by the rescattering of (tttc)j = j = o pions. Medium modifications 
are studied by coupling cr's to the nucleons by means of tadpole diagrams. Once inserted 
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into the model, the tadpole terms are found to contribute negligibly. Therefore, the only 
detectable nuclear medium modification on the (ini)i = j =Q system comes from the P-wave 
coupling of pions to p-h and A-h configurations. The model is able to predict the probability 
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Figure 8: Nuclear density distribution (p, full curve) as a function on the nuclear radius for 45 Sc. The 
probability of a tt2tt event to occur (H^n) is described by the distributions of the solid points: circles 320 
MeV (energy of incident pions), squares 280 MeV and stars 240 MeV. The 11^2^ distributions are normalized 
to unity, the U^tt and p yields are calculated with the model of Ref. [T5] . 

distribution of an interesting observable H^tt, which cannot be measured. Tl^-w is the com- 
posite probability of a ti2ii event to take place and be detected (by CHAOS) as a function 
of the nuclear radius. Fig. 8 shows the distributions for various energies of incident pions: 
circles 320 MeV, squares 280 MeV and stars 240 MeV. The distributions, when nor- 
malized to unity, vary slightly with energy thus indicating minor variations of the (average) 
nuclear density at which a n2n event occurs. In fact, p=0.36p n at T 7r +=320 MeV, p=0.37p„ 
at T 7r +=280 MeV, and p=0.35p„, at T^+=240 MeV. 



6. The T-dependence of the n + 



7T 



'tt^ 1 reactions in 45 Sc 



The experimental results are compared with the predictions of the model described in 
Sec. 5, and with the phase space of the 7r + 45 Sc — > tt + tt ± N 44 V reactions. Both model and 
phase space take the tt2tt CHAOS acceptance into account. The error bars reflect statistical 
uncertainties. Systematic uncertainties for T 7r +=280 MeV were discussed in Sec. 3. 



6.1 The 7T7T invariant mass, M. 



In order to check the consistency of the results, the distribution at 280 MeV is 
compared to the published M^ 1 cross section at 283 MeV [7j. The comparison is shown in Fig. 
9 top panels for the 7r + — > 7i + 7i ± reaction channels. The invariant mass of Sc (full diamonds) 
is normalised to the Ca invariant mass (open diamonds): as expected, the two distributions 
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display similar shapes since both kinetic energy and nucleus mass number are close. The 




Figure 9: Invariant mass distributions of pion pairs for the n + — > 7r + 7r ± reaction channels. The Ca data 
(open diamonds, [7J) were taken at 283 MeV and are used for normalisation. The Sc M TT data are indicated 
with full diamonds and were taken at 280 MeV. The phase space of the 45 S , c(tt+, ir+^N^X reactions 
are represented by the shaded diagrams. The curves are from |15l 1161 IT?] . The continuous lines denote the 
results of the full model. The dashed line in the upper left panel is obtained when pions are coupled only to 
p-h and A-h states. More details are given in the text. 

data (full squares) in the lower panels are the invariant mass distributions in deuterium, 
namely, of the elementary reactions 7c + n(p) — > 7i + n^p(p) and 7t + p(n) — ► 7T + 7i + n(n) 7 . The 
data are compared with the model predictions (full line) described in Sec. 5 Ref. (THj. The 
noticeable feature is that, for the reactions on 2 H, the calculations are able to reproduce 
the varying behaviour of the M n7T distributions at around threshold: M^^- is weakened to 
nearly zero while M n + n + reaches its maximum. 

In nuclei (Ca or Sc, upper panels) the situation changes. The 1^+^+ distributions resem- 
ble that of the elementary reaction and are well explained by the model predictions (solid 
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70 310 350 39D 430 

(MeV) 



Figure 10: Invariant mass distributions of pion pairs as a function of the incident pion kinetic energy. The 
phase space for 45 Sc(tt + , ir + TT ± N) 44 X reactions is indicated by the shaded diagrams. The full lines are the 
model calculations of |15l llfil IT?] . The kinetic energy of the pion beam is given in each panel. 



line). In the 7r + — > 7t + 7t~ channel, M nn displays a remarkable reshaping with respect to the 
elementary process, and the fully renormalised calculations (solid line) give only a partial 
account of the threshold enhancement of the data. Furthermore, the invariant mass intensity 
decreases slightly at around threshold when only the P-wave coupling of pions to p-h and 
A-h states is considered (dashed line). 

A further element of comparison is the reaction phase space (shaded diagrams in Fig. 9), 
which is normalised to data. Phase space describes the M 7r + 7r + distributions throughout the 
energy range, but largely misses the M„.+„- distributions. 

The features observed for Sc at 280 MeV are in general observed at all the examined 
energies. This is illustrated in Fig. 10, which shows the M n + n ± cross sections as a function 
of the incident pion kinetic energy. The tt + tt + invariant mass cross section is well explained 
by a microscopic model of the pion production reaction ^3], and the distributions are also 
described by the tt2tt phase space. This leads to the following conclusions. The dynamical 
traits of the (7r7r)7- =2 j = o interaction are understood, whether it occurs on the nucleon or in 
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nuclei. Furthermore, the nuclear medium does not appreciably change the irir interaction in 
the 1=2 J=0 channel with respect to the vacuum. 
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Figure 1 1 : Opening angle distributions of pion pairs as function of the incident pion kinetic energy (full 
diamonds). The phase space of the 45 Sc(ir + , 7r + 7r ± A'') 44 X reactions is indicated by the shaded diagrams. 
The full lines are the model calculations of ^JEHEI- ^he pi° n incident energy is given in each panel. 

In the 7r + — > 7T + 7r~ channel, the terms describing the irN — > ttttN amplitude must largely 
cancel out to explain the threshold depletion present in the model calculations. This view is 
also corroborated by the phase space behaviour, which approaches its maximum at around 
threshold. The enhancement measured at different kinetic energies (and nuclei) is only partly 
explained by theory, thus implying that either the cancellations are no longer effective in 
nuclear matter, or the a (i.e., strongly interacting (tttt) j = j =0 pairs) and its coupling to 
nucleons is not properly accounted for. In fact, other models of the (im) I= j =0 interaction 
in nuclear matter |UJ EI] find that the appearance of the a largely contributes to heighten 
the irir spectral function at threshold. These models, however, do not embed any pion 
production process. As a final comment, it is worthwhile noting that all the ir2ir data are 
well reproduced by both theory and phase space over the high invariant mass range. 



6.2 The irir opening angle, 
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For a reconstructed pion pair 7Ti7r 2 , the opening angle is defined as Q W17T2 = Q W1 - 0,r 2 , 
where the azimuthal angles G 7ri and Q n2 are measured in the lab system. The indices 1 
and 2 are assigned such that wi7r2 is restricted to vary from 0° to 180°. The tut opening 
angle is an observable that proved to be sensitive to the isospin state of pion pairs (see Sec. 
4 Fig. 7). The Q nn distributions are presented in Fig. 11 for all the examined energies, 
and compared with the reaction theory and phase space. In the 7r + — > 7r + 7r + channel, the 
Q nn distributions are generally flat over the angular interval 0°-180°. Such behaviour is 
predicted by the theory and phase space. The 0^+,,- distributions are peaked at around 
~ 0°, regardless of the energy. Such a feature is not described by the model calculations, 
which predict a smooth and decreasing intensity of Q nn when approaching zero degrees. 
For the sake of comparison, in Fig. 12, the Q nn distributions for the elementary reactions 



7r + n(p) — > 7T + 7T p(p) and ir + p(n) — > 7i + 7i + n(n) at T„ 
is characterised by a negligible strength below ~ 60 c 



=283 MeV [7] are shown. 0^+^- 
In this angular span, the strength 
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Figure 12: Opening angle distributions at 283 MeV kinetic energy. The shaded diagrams denote the phase 
space of the 2 H(tt + , tt + tt ± N)N reactions. The full lines are the model calculations of Refs. |15M16llT7| . 



detected in nuclei is a direct consequence of the medium modifications on the (7r7r)/ = j =0 
interaction, since the ttN — > tcttN reaction is a quasifree process. The tt + — > -k + tt + channel 
behaves differently. The elementary reaction already exhibits strength below 60°, which 
does not increase in Sc. Theory [TH] is able to predict both the intensity and the shape 
of the 0^+^+ distributions for Sc. For 2 H, calculations account correctly for the measured 
intensities over the interval 0° < 0^+^+ < 180°, but the distribution is predicted to peak at 
about 60°, 20°-25° below the measured value. 



For the tttt opening angle analysis, the conclusions that can be drawn are similar to 
those for the invariant mass. At intermediate energies, nuclear matter clearly modifies the 
(tt7t)j = j =0 interaction, while it leaves the (7r7r)/ = 2j=o interaction rather unaltered. The theory 
is able to explain the n2n process only when it occurs in the 1=2 J=0 channel. 



6.3 The 7r kinetic energy, T 
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The observable under consideration, and illustrated in Fig. 13, is the kinetic energy of 
a single pion. In the 7r + — > 7r + 7r~ channel, positive (full diamonds) and negative (open 




Figure 13: Single pion kinetic energy distributions as a function of the incident pion kinetic energy. Positive 
pion distributions are indicated with full diamonds, negative pions with open diamonds. The shaded diagrams 
are the result of phase space simulations for the tt + 45 Sc — > n + n ± N 44 X reactions. The full lines are the 
positive pion energy distributions which are derived from the model of Refs. [1511161 IT?] . 

diamonds) pion distributions can barely be distinguished throughout the energy range ex- 
amined. The full lines [TS] describe the positive pion energy distributions. Calculations 
provide a fair description of the measured distributions, but tend to overestimate the inten- 
sities of the data in the high-energy range. They predict greater abundances of energetic 
pions, which are also in excess of those expected from 7r + 45 Sc — > tt + tt~p 44 Sc phase space 
(shaded diagram). The T 7r + distributions are well accounted for by the 7r + 45 Sc — > n + n + n 
A4 Ca phase space simulations, which are normalise to the n2n data. On the theoretical side, 
calculations for the 7r + — > 7r + 7r + channel provide a degree of agreement with data similar to 
that found for the 7r + — > 7T + 7r~ channel. 



The study of the T w observable completes the present discussion. It confirms the frame- 
work already sketched out by the study of the previous observables, i.e. and 8„. The 
7r + — ► 7r + 7r + reaction channel can be explained in terms of both model calculation and phase 
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Table 2: Experimental total cross sections of the tt + — > tt + 7t~ (H — ) and 7r + — > 7r + 7r + (++) 
reaction channels on 45 Sc as a function of the pion incident energy T n +. The uncertainties 
associated to the phase — space cross sections are omitted since they do not exceed the 
uncertainties of the linear cross sections. In the last two columns the a^" /(?t + ratios are 
reported for Sc and H, respectively 



\MeV] 




a T [lib] 








7<4 + 




linear 


phase space 








+- 


++ 


+- 


++ 


45 Sc 




240 


1225±164 


143±19 


1209 


141 


8.6±1.6 


7.1±1.2[25J 


260 


2595±348 


347±47 


2559 


341 


7.5±1.4 


6.1±0.9j23 


280 


4263±571 


507±68 


4233 


505 


8.4±1.6 


7.4±1.0j2S] 


300 


6660±892 


640±86 


6643 


637 


10.4±2.0 


9.6±1.5|2Ej 


320 


9160±1227 


826±110 


9148 


824 


11.1±2.1 


12.4±3.0 26J 



space. The n + — > 7r + 7r~ data are only partially clarified; in the present case, the model 
predictions fail to describe the high-energy tail of the T w distributions, which follow phase 
space. 

6.4 The ir2ir total cross sections, <jt 

With the CHAOS spectrometer one can directly determine d?a / ' dO^dVL^dVt^, where O n7T 
represents a generic observable, dfl v is the solid angle into which a charged pion is scattered 
and dQ n = cfcos(0 7r )<i ( lv In order to assess ax, the multiple differential cross section has 
to be summed over and Q^s. The limited ^-acceptance of CHAOS, however, allows 
the determination of da/d®^ only at = 0° and = 180°, where = = 
$7r 2 |" Within this interval, the $— dependence of the h2tt reaction is presently assumed 
either to behave linearly, or to follow the 45 5'c(7r + , n+ix^NY^X phase space (see Sec. 3, and 
Fig. 4 for details). In this framework, the results obtained for a? ar e those reported in Table 
2. In the last two columns, the ratios /a^ + for the two reaction channels are listed. For 
a given T^, the ratios for Sc and H agree within the error bars, which implies that a similar 
reaction mechanism underlies the 7r27r reaction whether it occurs in the nucleon or in nuclei. 
This confirms that the ti2h reaction in nuclei is a quasifree process. 

7. Existing n2ix results 

Two collaborations have recently presented their studies of tttt dynamics in nuclear mat- 
ter in proximity of the 2m 7r threshold. These are the Crystal Ball (CB) Collaboration at 
the AGS and the TAPS collaboration at MAMI. Comparisons can only be made with the 
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general features of the data, since pion pairs are produced by different reactions, or different 
isospin channels are studied, or different energies and nuclei are examined in the different 
experiments. 

The CB Collaboration has presented results on the tt~A — > 7i°7i°A' reaction for A: H, 
12 C, 27 Al and 64 Cu, at an incident pion energy of T n - =291.6 [22] • Neutral pion pairs have 
1=0, therefore they can be directly compared to 7r + 7r _ pairs A recent article 

j2H] compares the CB and CHAOS data, pointing out the relevant common features shared 
by the two data sets. The observable is of particular interest in the comparison. 
is proportional to the ratio of to [H] d ■ Consideration of such a ratio lessens the 
importance of the difference in the CB (~93% of 4tt) and CHAOS (~12% of 4ir) acceptances. 
Comparison shows that the behaviour of C^ 0jT o and C^ +7r _ is similar over the entire energy 
range. 

mi dynamics in nuclear matter has been also studied by the TAPS Collaboration. TAPS 
has used the photoproduction reaction 7 — > nn on H, 12 C and 208 Pb, at photon energies 
from 400 to 460 MeVj21]. At these energies, gammas can penetrate to the interior of nuclei, 
therefore probing nuclear densities (p ~ 0.35p n for 12 C and p ~ 0.65p n for 208 Pb) higher than 
the pions in the present work (p ~ 0.36p n ). In the case of TAPS, medium effects on the 1=0 
7i7x interacting system are understood to be broader than the CHAOS ones, but the nuclear 
distortions on the two pions being detected are also broader, which has an impact on the 
cross section. Therefore, a simple comparison of C^ 07r o to C^ +7r _ may be misleading, since 
distortions do not affect M^.. TAPS has redefined the composite ratio C nn = M™/M% n . 
With this new definition, the TAPS C nn agrees well with the CHAOS C wn (Ref. [21]). In 
addition, strong medium modifications have been observed only for the 7 — > 7i°7i° reaction 
channel; that is, for pion pairs interacting in the I=J=0 channel. Medium effects on the 
isospin 1=1 7i7i pairs have not been detected. Further points of contact between the TAPS 
and the CHAOS results are discussed in Ref. [21]. 

8. Conclusions 

This article presents the results of an exclusive measurement of the pion-production 
7r + — > 7r + 7T ± reactions on 45 Sc, at incident pion energies of 240, 260, 280, 300 and 320 MeV. 
In previous measurements the same reaction has been examined on 2 H, 4 He, 12 C, 16 0, m Ca 
and 208 Pb nuclei, at an incident pion energy of ~280 MeV [UI21EII11E1EIIZ]- A two-fold 
interest has been driving these measurements: understanding the general traits of the pion- 
production reaction in nuclei, and studying the tttt dynamics in nuclear matter. To this end, 
the CHAOS magnetic spectrometer at TRIUMF has been used, which permitted the analysis 
of 7r + 7r + and 7r + 7r~ pairs simultaneously. The study of the medium modifications of the tt + 7t~ 
interacting system is of central importance, being the lightest system which carries the 1=0 
J=0 quantum numbers, that is, the quantum numbers of the QCD vacuum. At the present 



19 



energies, the 1=2 pion pairs interact very weakly, thus, only weak medium modifications on 
the 7r + 7r + system are expected. Therefore, the comparison of the two isospin channels yields 
direct indications of the nuclear medium effects on the (7T7r)/ = j = o system, that is, on the a 
meson. 



The A- and T-dependence studies of the 7r27r reaction lead to the following conclusions. 



1. The 7t2it reaction in nuclei is a quasifree process, which involves a single nucleon 
ttN — > ititN. 

2. An analysis of the p n7T — dependence of reveals that near the 2m 7r threshold tttt 
pairs carry the highest momenta available, and the opening angle distribution of -k + tt~ 
pairs is peaked at ~0°. This behavior is not explained by the reaction phase space or 
by theory. 

3. The 7r27r data are compared with a complete model of the pion-production reaction[15j. 
At the present energies, calculations indicate that the average nuclear density at which 
the process occurs is p ~ 0.36p n , that is, the process is confined to the nuclear skin. 

4. The observables M nn , Q n7T and T n have been carefully examined and compared to 
model calculations and reaction phase space, both of which take the CHAOS accep- 
tance into account. The dynamics of the 7r + — > 7r + 7r + channel is well explained by 
theory. Phase space provides a good description of the kinematic behaviour of the re- 
action over a wide range of energies. The nuclear medium does not detectably modify 
the (nn) t=2 j=o interaction. For example, the threshold behavior of M 7r + 7r + depends 
negligibly on A and T. This also indicates that pion final state interactions cannot 
account for the enhancement observed at around the 2m,,- threshold. The dynamics of 
the 7r + — > 7r + 7r~ channel is far more varied. For the elementary ttN — > ttttN reaction, 
the M n + n - strength at around threshold is depleted over a wide range of incoming pion 
energies, from 223 to 305 MeV Ref. [23] and this work Fig. 9. Thus the terms forming 
the elementary amplitude must cancel out. In nuclear matter, M n + n - rapidly gains 
strength, which is not caused by the final state interactions of pions with the resid- 
ual nucleus. Therefore, either the medium removes cancellations or it produces new 
events. According to the authors of Refs. EH], an event is the appearance of the a 
meson, a genuine QCD system formed by a qq state, which largely reshapes its spectral 
function in response to nuclear matter. The M 7T + n - threshold strength is primarily due 
to this occurrence. However, this basic result is not corroborated by a model of the 
7r27r reaction. The model used in this article accounts for the pion-production process 
[To] . At around threshold, M^^- builds up strength via the P-wave coupling of pions 
to p-h and A-h nuclear excitations. This process also underlays the dynamical origin 
of the a meson: a tttt interacting state in the I=J=0 channel. This resonant state is 
found to drop in mass and width at finite densities, i.e. at p < p n [27|. However, the a 
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strength at around threshold is insufficient to explain the data. This specific issue de- 
serves the attention of theorists, since it is of primary importance in the understanding 
of whether the invariant mass shift toward threshold is due to the partial restoration 
of chiral symmetry or to other essential events. Standard nuclear effects are already 
accounted for by the model used in the present analysis. This general conclusion also 
applies to the Q nn and T^ observables. 

5. The 7r27r results obtained by the CHAOS Collaboration have been compared to the 
results of the CB and the TAPS Collaborations, which are the only available ones. 
The results discussed in the CHAOS articles share relevant common features with the 
results presented by these two collaborations. The discussion over the composite ratio 
= M^/M^ has been only briefly developed since it is the subject of a forthcoming 
article. 
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